Evolution of cosmological event horizons in anisotropic universes 
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We study the evolution of cosmological event horizons in anisotropic Kasner universes in the pres- 
ence of a positive cosmological constant by analyzing null geodesies. At later times, the asymptotic 
form of cosmological horizons is the same spherical surface as the de Sitter horizon. At the early 
times, however, it has non-spherical shape with its eccentricity decreases with time. The horizon 
area increases with time respecting the second law of thermodynamics. The initial shape of the 
cosmological horizon takes the form of a needle or pancake surface depending on the nature of 
the background spacetimes. We also discuss that the presence of the holographic dark energy will 
modify significantly the initial evolution of the anisotropic universes. 
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Following the Copernican principle [l[ , the physical properties of the universe are spatially homogeneous viewed on 
. a sufficiently large scale. This is consistent with nowadays experiments up to of order 10 -5 . Based on the principle, 
■ the Robertson- Walker spacetimc with inflation at early stage, on which small perturbations evolve to give cosmological 
r ^ , structures, is widely used to describe the present universe. However, there is no reason to believe that our universe was 
1™^ ■ initially homogenous and isotropic. Gcncrically, it could be highly inhomogeneous and anisotropic. The cosmic no-hair 
t — \ theorem Q shows that in the presence of inflation most of the initially anisotropic cosmological spacetimes become 
isotropic rapidly if the matter fields residing in the universe satisfy the dominant energy condition. The theorem 
1 also implies that it is inevitable to have a highly anisotropic pre-inflationary universe if any small anisotropy remains 
q-i after inflation. Since the anisotropy will be exponentially enhanced if we observe the universe backward in time 
^ " and we cannot make that process stop by adding ordinary matters or radiations. The vestige of the pre-inflationary 
j^Tj anisotropics can be seen through the large scale anomalies in the WMAP power spectra if the e-folding of inflation 
'— is not larger than 64 for GUT scale inflation. [H, H|. In fact, there is possibility that the effect of pre-inflationary 
anisotropics has been already detected in the Cosmic Microwave Background (CMB) temperature anisotropics, as an 
apparent alignment of the CMB multipoles on very large scales, so-called "Axis of Evil" and giant rings [Tlj . 
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The standard cosmological model is the A-CDM model where most part of the universe is composed of the cosmologi- 
OO ■ cal constant and cold dark matters. The spacetime has a future event horizon located at a distance Rh — a(r) 



spy' 

where a(r) is the scale factor of the Robertson- Walker metric and r is the comoving time. In the presence of asymp- 
totically accelerating expansion of the universe, any observers will be surrounded by the future event horizon, which 
restricts the observers access to the information of the universe. This lack of information is represented by a kind of 
entropy proportional to the area of the cosmological horizon divided by the Planck area similar to the black hole en- 
tropy. The generalized second law of thermodynamics for the Robertson- Walker spacetime was proved by Daives [l2[ 
and Pollock and Singh [l3| showing that the total entropy of gravity+matter does not decrease through physical 
processes. It would be interesting to find the cosmological horizon of anisotropic universes whose geometry are not 
isotropic and investigate the possibility to generalize the second law of thermodynamics to such spacetimes. At the 
j | 1 present universe, the cosmological horizon is deeply related to dark energy which dominates the later accelerating 
d \ expansion of the universe. It is well known that the portion of the dark energy is increased during the matter dom- 
inated regime and is negligible at the radiation dominated regime. This is why one usually neglect the dark energy 
contribution during inflationary and pre-inflationary regimes of the universe. However, in the case of an holographic 
dark energy, one cannot neglect the dark energy portion at the pre-inflationary regime. If we observe the universe 
backward in time during the inflation, we may see that the portion of the holographic dark energy density will be 
exponential increased until it is comparable to the inflaton energy density at the pre-inflationary universe as shown 
m Refs. [3 [la ]. Therefore, it is reasonable to introduce the holographic dark energy on the anisotropic inflation 
models and ask the possibility that the behavior of anisotropy in Kasner models is modified by its presence. In fact, 
this gives a possibility to cure the pathological initial condition problem of the gravitational wave perturbation fl(| , 
which happens because of the initial anisotropic singularity of Kasner models. 

We are interested in the homogeneous anisotropic cosmological solutions of Einstein gravity with a positive cos- 
mological constant, A(= 3i? 2 ), in the absence of matter fields. The spacetime solutions are known as the Bianchi 
universes, which are divided into nine classes. In this paper, we study the simplest case: the Kasner spacetimes with 
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two dimensional plane symmetry with metric, 

ds 2 = -dr 2 + a 2 (r) (dx 2 + dy 2 ) + c 2 (r)dz 2 , (1) 

where a(r) and c(r) arc independent scale factors of the planes and of the direction orthogonal to the plane, respec- 
tively. The space-time model is a good testing ground in analyzing the properties of anisotropic universes since it is 
simple and has various important features of the whole anisotropic universes. 

There are two different solutions of distinct features. One is given by ds 2 := ds 2 with scale factors 

a(r) := a r (r) = a sum* (3.HY) tanh~' > C ( T ) ;= Cl '( T ) = c o sinh ' ( SHt ) tanh^ (— jj— ) • ( 2 ) 

The spacctimc has a Rindlcr-like event horizon at r = since a r approaches to a finite value and c r goes to zero 
linearly. The spacetime is connected to the region with negative r through the horizon where a timelike singularity 
is located. As shown in Ref. [l7|, the geometry is closely connected with the anti-dc Sitter black brane solution 
discovered by Lemos [l8|. The other solution, ds 2 := ds 2 , has an initial singularity and their scale factors are given 
by 

a(r) := a s (r) = a sinh^ (3ifr) tanh^ (— ) ; c ( r ) : = c s(t) = c sinh^ (3-Hr) tanh~3 (j^-^j ■ (3) 

The metric component a s goes to zero and c s diverges as r — > 0. As a result, there appears a spacelike singularity, 
which is present in most Bianchi spacetimes other than ds 2 style. The solutions behave initially as Kasner spacetimes 
in the absence of cosmological constant followed by a de Sitter phase, as discussed in In this paper, we set 

ao = 1 = Co, which can always be done by an appropriate rescaling of the spatial coordinates x l (i = 1,2,3). 

Now we consider the null geodesies on this spacetime. Consider a flash of light with conserved comoving momentum 
Pi '■= |p| (sin cos sin sin 0, cos 0) which departed from the origin at time r, where 9 denotes the angle of the 
comoving momentum with respect to the z direction. At time To > r, the light will arrive at the point, 

r l = (x, y,z); x(t, to) = sinfl cos f(6>, r , r), y(r, r ) = sin6> sin0 f((9, r , r), z(t, t ) = cos6» g(0, r , r) (4) 

where 

r To c(t') r To a(r') 

mro,r):= d T ' 1 ) q(6,t ,t) := d T ' a[ > .(5) 

Jt a(r')y sin 2 9 c 2 (r') + cos 2 6 a 2 (r') Jt c(t')J sin 2 6 c 2 (r') + cos 2 9 a 2 (r') 

The angle 9 docs not directly related to the physical zenith angle of the points (x, y, z) but is related to the conserved 
momenta. 

Note that, for large r, the two scale factors will eventually evolve exponentially and become the same as that of 
the de Sitter spacetime. In this sense, the existence of the cosmological event horizon is guaranteed. As tq — > oo, the 
light will arrive at the spacial coordinate 

ri{T,9,(j)) = lim r l (T,T ). 

T ->OO 

Conversely, a flash of light with comoving spatial momentum — pi which departed from the event (r,Xh,yh, Zh) will 
arrive at the origin as tq — > oo. This clearly implies that the surface, 

Sj{r) ee {(r, x h , y h , z h ), < 9 < tt, < < 2tt, t > 0}, (6) 

forms the cosmological event horizon of the universe at the comoving time r. If a = c, the formula reproduce the 
distance to the cosmological event horizon in the Robertson- Walker spacetime. 
Apparently, the horizon surface satisfies 

x 2 + y 2 z 2 



f 2 (6>, oo,t) 2 (6>,oo,t 

with f and q dependent on time and angle 9. If the two functions f and g arc independent of the angle, the surface 
forms an ellipse. This case happens only asymptotically as r — ¥ oo. Explicitly, the distance to the cosmological event 
horizon is dependent on the polar angle 9 and is given by 



R h (r,6) = JR 2 + R 2 ± , (7) 
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where the distance along the parallel and orthogonal directions to the x-y plane are given by 

R\\(t,0) = a{T)^x 2 h {T)+yl{T) = sin0a(T)f(0,oo,T), R±{t, 9) = cos 9c{t)q{9, co, t). 

In this paper, we explicitly calculate the distances to cosmological event horizons for the two metrics (0) and (J3J. Let 
us first study the initially regular spacetime, ds 2 , case. Considering the perturbation around the Kasner spacctimc 
with positive cosmological constant, the gravitational wave modes rapidly grow and have instability at the initial 
Kasner regime, gencrically. Only for the present solution, the instability can be suppressed to give a well defined 
initial condition for the gravitational perturbations JiJ. Eq. ([5]) is intcgrablc to give a closed form and the spatial 
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FIG. 1: Shape of the event horizon (HR^,HR T ± ). Here, 3Ht/2 — 2.5, 1.5, .5, .1, and .001, respectively from small eccentricity. 
We choose H = 2/3. 



distance to the event horizon for a given angle 9 becomes 
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g x (sin 9, cosh (^^-j] 



1/2 



(8) 



where 



f r ( X ) = 1^/3 / dx > 



2 2/3 
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v^r(i/3) 



- y/x 2 -l 2 F X (1/2, 5/6, 3/2,1 



2r(5/6) 

is a monotonically decreasing function of x (> 1) with / r (l) = ^Y{b/&) an< ^ um f?( x ) = 1- The function, 

, (2/ + l) 1/6 
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monotonically increases from zero to one as y (> 1) increases. 

For the limit r — > oo, the distance to the horizon approaches to that of the dc Sitter space as R r h — > 1/H, and the 
cosmological horizon forms a spherical surface. On the other hand, around the time r — > 0, the z-direction distance 
to the horizon goes to zero as 



rlog 



3Hi 



Note that this distance is independent of the angle 9. As r — > 0, the distance to the horizon parallel to the plane 
approaches to a finite number, 



Rn 



0F| sin 1 / 3 ^ T(l/3) 
3H T(5/6)' 



Therefore, the horizon take the form of a pancake surface of infinitesimal thickness as shown in the case of 3Ht/2 
0.001 in Fig. [I] The area of the initial event horizon, 



r \2 



A = 2n(m) 



2tt 2 T 2 (l/3) _ 12.3535 
9lPr 2 (5/6) ~ H 2 ' 



is slightly smaller than that of the final dc Sitter horizon Att/H 2 ~ 12.5664/i/ 2 . The horizon area are increased 
through the evolution. 

We next consider the spacetimcs with initial spacelikc singularity. Even though, the metric ds 2 is a specific example 
of general Kasncr spacetimcs, it still has most important features of the initially singular anisotropic spacetimcs. With 
metric (0, Eq. §5§ is integrable to give a closed form and the spatial distance to the event horizon for a given angle 
9 becomes 



H 



fs( 



sinh(2f^ 
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<? s (cos#, sinh' 
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(10) 



where 



/.(*) = 



2x 2 / 3 



dx' _ 1 1 4 1 
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monotonically increases from zero to one as x increases from zero to infinity. The function 



3 (c, y) 



(y + i) 1 ' 2 f 00 , , y' 1/e A 
32/ 1 / 6 J v y (y' + l) 3 / 2 \ y' + l, 

(y + H /2 V ^(n + l/2) 2n ( _ )5/3+ ^ (1/3 + _ 1/2 _ n) 
3y i/6 Z_^r(l/2)r(n + l) 7 y ' 



-1/2 



(12) 



monotonically decreases from infinity to one as y (> 0) increases. As a result, the distance to the horizon asymptotically 
approaches to that of the de Sitter space as R s h 1/H. 

On the other hand, around the initial singularity, the z-direction distance to the horizon diverges as 



r(4/3)r(l/6) cos 6*2^(1/3, 1/2, 3/2, cos 2 9) 



3sinh l/3(Mr ; 



The distance to the horizon parallel to the plane approaches to zero as 

r(4/3)r(i/6) s . nl/3gs . nh2/ 3 /affr\ 

11 V^H \ 2 J 

Therefore, the initial horizon take the form of a needle surface of infinitesimal thickness and infinite length as shown 
in Fig. [2j The area of the initial cosmological horizon goes to zero as 



A - 27ri?fji?l oc 



sinh 1 / 3 (Mi) 
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FIG. 2: Shape of the event horizon (HRl, HKj_). Here, 3Ht/2 = 1.5, 1, .5, .1 and .01, respectively from small eccentricity. 



The horizon area is increased through the evolution. We plot the event horizons at several cosmological times in 
Fig. [2] The whole form of the event horizon can be obtained by rotating the graph around the z-axis (the horizontal 
axis). For Ht 3> 1 the event horizon becomes a spherical surface. The eccentricity is very high initially even though 
its shape is not an exact ellipsoid. 

Let us summarize the results in this paper. At the future infinity, the cosmological event horizon has spherical 
shape and is located at 1/H. For times other than the future infinity, it is non-spherical and has non-zero eccentricity, 
which is maximized at the initial time. In the case of the initially singular Kasner spacetimes, it was shown that the 
initial cosmological horizon takes the form of a needle surface with infinite length and infinitesimal thickness. The 
initial surface area goes to zero since the thickness decreases faster than the increase of the length if we observe the 
universe backward in time around the initial singularity. In the case of the initially regular Kasner spacetime, the 
initial cosmological horizon takes the form of a pancake surface on the symmetric plane with infinitesimal thickness. 
For both cases, the surface areas of the cosmological event horizons increase with time. Therefore, the 2 nd law of 
thermodynamics holds in the anisotropic universe allowing the entropy interpretation of the horizon area. 

Calculating surface gravity on the horizon is a challengeable task since the Killing vector generating the horizon is 
not known. The vector following the horizon for a given value of 9 and <fi, are given by 

a dxf t _ c sin 9 cos 4> csinfl sva<f> acos9 . 

^ r a \J cos 2 9 a 2 + sin 2 9 c 2 a\J cos 2 9 a 2 + sin 2 9 c 2 c\J cos 2 9 a 2 + sin 2 9 c 2 

One can easily show that this is null v a v a = but not a Killing vector. The horizon generating Killing vector \a will 
be proportional to v a at the horizon. We defer this subject of finding the vector to later work. Rather, we roughly 
estimate the surface gravity near r ~ by using analogy to the de Sitter spacetime. Near the cosmological horizon 
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located at x = H x e Ht with y = = z, the metric of the de Sitter spacetime becomes 

ds 2 ~ -2a(T)dTd£ + a 2 (r)(dC 2 + dy 2 + dz 2 ), 

where £ = x — H~ 1 a~ 1 and a(r) = e Hr . The surface £ = define the cosmological horizon. The surface gravity in 
this spacetime is given by the derivative of the logarithm of scale factor, a -1 (da / 'dr) . In this calculation, the direction 
y and z are irrelevant. 

In the anisotropic Kasner spacctimes, we may proceed in a similar manner to estimate the surface gravity of the 
cosmological horizons at early times. For the case of the initially singular metric ds 2 , most part of the horizon is almost 
parallel to the z-axis initially since the horizon has needle-like shape. Therefore, we may assume rough translational 
symmetry along z. We estimate the surface gravity at the event (t <C H , x^, Uh = 0, Zh ~ 0) on the horizon. Since 
the geometry is almost needle-like, we may simply calculate the surface gravity at 9 = tt/2. Similar calculation as 
above leads to k s ~ H a = a/a oc 1/r at the early times. For the case of the initially regular metric ds 2 , the horizon 
is almost pancake shape located on the x-y plane initially. We may estimate the surface gravity at the point 6 = 0. 
Similar calculation leads to /c r ~ H c = c/c oc 1/r. 

The surface gravity does the role of a Hawking temperature for the geometry. Usually, a kind of thermal energy can 
be related to the horizon and at the present case one usually relate it to the holographic dark energy 0, Hl| which 
explains the cosmic expansion of our universe. It was also pointed out that the second law of thermodynamics holds, 
dA J^ H > 0, for the case of a Robertson- Walker cosmology with holographic dark energy. In addition, the second law 
can be used to determine when the dark energy dominated era begins in the universe |22j . 

Let us assume that there present small amount of holographic dark energy and matter in the later (but before 
inflation) evolution on the background of Kasner spacetime. Then, let us see how does the energy densities behave 
as we go backward in time. The holographic dark energy density at earlier times behaves as 

T H S kA 

Ph 



V in R 2 (t,it/2)R ± (t,0) : 



where Tjj, S, and Vin denote the Hawking temperature, entropy, and volume of the universe inside the horizon, 
respectively. Therefore, the holographic dark energy density in the initially regular or singular spacetimes become 

H 2 / 2 V 2 H 2 
P R S ~ iTj ^ lo s Try ' p H , s 



(Ht) 3 \ "SHtJ ' ' "' s (Ht) 7 / 3 ' 

The energy density of matter field behaves as p m ^ ~ po / (a 2 (r)c(r)) ~ p /(Hr) and p m .s ~ Pol '(Ht) 2 / 3 , respectively. 
As time goes back, the holographic dark energy density increases much faster than that of the matter. Therefore, 
the holographic dark energy eventually dominates the universe at some early times. If one consider the reaction of 
the holographic dark energy to background geometry, which is not taken care of in this paper, the evolution of the 
Kasner universe should be modified. 
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